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(57) Abstract: A method of processing 
marine magnetic data is disclosed which 
comprises towing first and second 
magnetometers Mx and M2 behind a 
ship (10). Raw magnetic gradient data 
is obtained from the sensors and the 
trend of the gradient of the ship bias 
detected by the sensors determined. The 
method includes subtracting the trend 
from the raw magnetic gradient data 
to obtain corrected gradient data, and 
processing the corrected gradient data 
to provide a data output. 
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METHOD OF PROCESSING MARINE MAGNETIC GRADIENT DATA AND 
EXPLORATION METHODS USING THAT DATA 



Field of the Invention 
5 This invention relates to a method of processing marine 
magnetic gradient data and also to exploration methods 
using that data. 



Background Art 

10 Exploration of the earth to locate oil, minerals or other 
valuable deposits, uses seismic techniques in order to 
provide information about the subsurface structure of the 
earth so that an assessment can be made as to whether a 
particular deposit may or may not be present. In 

15 particular, seismic techniques provide an indication of 
various subsurface structures including porous sandstone 
or fractured carbonates which may contain oil or other 
valuable deposits. 

20 Whilst this technique does provide an indication of 
subsurface structures, it is possible a particular 
structure which may appear to be of interest, in fact 
represents a significant hazard to drilling. For example, 
it is possible that subsurface volcanos exhibit very 

25 similar characteristics to anticlines which may contain 

oil. If a drilling operation is initiated and it is found 
that the anticline is in fact a volcano, then a 
considerable financial loss results because of the cost of 
establishing the drilling platform. This is particularly 

3 0 the case in relation to marine exploration because of the 
significantly higher drilling cost. 



35 



In order to provide more information concerning the nature 
of subsurface structures, magnetic data of the exploration 
area is obtained. 
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A magnetic gradient survey enables magnetic data relating 
to the survey region to be obtained and this magnetic data 
can be used to provide information relating to the nature 
of subsurface structures. If the magnetic survey is 
5 overlayed with the seismic survey, then structures which 
appear to be of interest from the point of view of the 
seismic survey, can be further considered in the light of 
the magnetic data so that a clearer indication can be 
formed as to whether the structure is a structure which 

10 may contain a deposit of interest such as an oil deposit, 
or whether the structure exhibits magnetic phenomena, such 
would be the case with a volcano. Thus, the location of 
drilling platforms can be decided with more precision to 
avoid subsurface structures which could be hazardous from 

15 the point of view of a drilling operation. However, 

conventional techniques for processing magnetic gradient 
data contain considerable distortion and anomalies due to 
unwanted magnetic effects including the ship bias referred 
to above. 

20 

The conventional method of obtaining that data in relation 
to onshore exploration is to tow magnetometers behind an 
aircraft to obtain magnetic data in relation to the survey 
area. In the case of marine exploration, the 
25 magnetometers are towed behind a ship. 

The usefulness of marine magnetic data which is obtained 
in this manner has been limited by the quality of the data 
which is obtained. One of the major problems with marine 
30 magnetic data acquisition is interference or so-called 

ship bias which is created by the magnetic field induced 
by the ship which tows the magnetometers. 

In conventional marine magnetic data acquisition, two 
35 magnetic field sensors, which are generally referred to as 
fish, are towed behind a ship. The magnetometers are 
connected to a tow line and the magnetometer which is 
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closest to the ship is towed at a distance of some 300-600 
metres behind the ship to avoid the effect of ship induced 
magnetic field. Furthermore, the magnetometers are 
separated by a distance of in excess of 100 metres* The 
5 reason for the length of the tow line and the separation 
of the magnetometers is to reduce ship bias and therefore 
provide data which is relatively free of that bias. 



However, because the towing line is so long and the 
10 distance of separation so great, the degree of drift of 
the magnetometers in the ocean as the magnetometers are 
towed behind the ship is considerable. Furthermore, the 
processing techniques used to acquire magnetic gradient 
data assume that after a given time period, the trailing 
15 magnetometer will be towed to a position which coincides 
with a previous position of the magnetometer which is 
closest to the ship. This assumption is made during the 
processing of the magnetic data. However, because the tow 
line is extremely long and the distance of separation 
20 between the magnetometers is considerable, drift of the 

magnetometers due to sea currents and the like means it is 
unlikely that the trailing magnetometer will actually 
occupy the same position as the first magnetometer after 
that given time period. 

25 

Furthermore, in the acquisition of the magnetic data, the 
ship is required to travel along predetermined survey 
lines and because of the drifting of the sensors, the 
sensors do not actually travel along the survey line. In 
3 0 the southern hemisphere, if the magnetometers drift to the 
north of the survey line, the recorded gradient data 
between the two magnetometers will have larger distortion 
than if drifting towards the south of the survey line. 



35 



Summary of the Invention 

The object of the invention is to provide an improved 
processing method which accounts for ship-induced bias and 
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which therefore enables the magnetometers to be towed by a 
tow line which is much shorter than the conventional tow 
line and also be located closer together to reduce the 
amount of drift and also increase the likelihood that the 
5 trailing magnetometer will occupy the same position which 
was previously occupied by the leading magnetometer after 
a given time period. 

The invention may be said to reside in a method of 
10 processing marine magnetic data obtained by towing first 
and second spaced apart sensors behind a shipT the method 
including: 

obtaining raw magnetic gradient data from the 

sensors ; 

15 determining the trend of the gradient of the ship 

bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
2 0 a data output. 

Because the method determines the trend of the gradient of 
the ship bias, a more accurate estimation of the actual 
gradient bias introduced by the ship, having regard to any 

2 5 drifting of the sensors, is taken into account. By 

subtracting the trend of the gradient bias from the raw 
gradient data, the bias introduced by the ship can be 
removed from the data so that more accurate data is 
obtained. Because the present method eliminates the 

3 0 effect of the ship bias, the sensors can be towed much 

closer to the ship and at much closer separation 
distances, thereby reducing the amount of drift of the 
sensors and therefore resulting in the sensors being more 
likely to follow the actual survey line upon which the 
3 5 ship travels and also for the rear sensor to occupy the 
same position at a later time which was previously 
occupied by the front sensor. Thus, the equality of the 
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output: data is greatly improved and therefore provides 
data which gives a much more accurate indication of the 
subsurface magnetic characteristics of the survey area, 
which can be then used in combination with seismic data to 
5 make an assessment as to the viability of further 
exploration or a drilling, program. 

Preferably the method includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
10 magnetic gradient data obtained by the sensors, and 

determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

Preferably the raw magnetic gradient data is calculated as 
15 the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors. 

The estimate of the gradient of the ship bias may in fact 
20 simply be derived from the raw magnetic gradient data 

which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

25 Preferably the trend of the gradient of the ship bias is 
determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

Preferably the processing comprises the application of the 
3 0 corrected gradient data over sample intervals, integration 
of the corrected gradient data into total magnetic 
intensity data, and application of a low pass filter to 
the integrated total magnetic intensity data to provide 
the data output. 

35 

Preferably the total magnetic intensity is obtained by 
integrating the corrected gradient data in accordance with 
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the following equation: 

TML= £ [G c Ax/(0] + Ms*, 

wherein 6 C is the corrected gradient data 
5 obtained after subtraction of the trend of the gradient of 
the ship bias from the raw gradient data, and 

M fl tat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
10 other. 

The values, of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter* 

15 Preferably the gradient of the ship bias is calculated 
according to the following equation: 

M f {x)-M r {x) _ [D(r, ) - D(r 2 )] + [M b (t } )-M h (t 2 )] 
Gbias = II = Al 

2 0 where : 

Af/jc) = M e {x) + D(t } ) + M b {t x ) 
M r (x) = M e (jc) + D(r 2 ) + M 6 (r 2 ) 

25 

where M £ is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Mo(x), diurnal variation D(t x ), ship bias M b (ti) caused by 
ship-induced field, sensor drifting and sensor 

3 0 perturbation, all at time t x and along track distance x, 

and at some later time t 2 # Al is the distance between the 
sensors, D(t a ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
35 rear sensor M r , 
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In one embodiment of the invention the first and second 
sensors which are towed behind the ship are included in a 
group of three or more towed sensors. 

5 In this embodiment, preferably the number of sensors towed 
behind the ship comprises three sensors. 

According to this embodiment data from any two of the 
sensors is used to provide the raw magnetic gradient data. 

10 

However, the embodiment may obtain raw gradient data from 
all three sensors and determine the trend of the gradient 
of the ship bias detected by all three sensors. 

15 The invention may also be said to reside in a method of 

obtaining gradient data for an exploration program, the 

method including: 

towing first and second sensors behind a ship 

along predetermined survey lines; 
20 obtaining raw magnetic gradient data from the 

sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 

subtracting the trend from the raw magnetic 
25 gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
a data output. 

Preferably the method includes the step of determining an 
30 estimate of the gradient of the ship bias from the raw 
magnetic gradient data obtained by the sensors, and 
determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

35 Preferably the raw magnetic gradient data is calculated as 
the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
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divided by the distance between the sensors. 

The estimate of the gradient of the ship bias may in fact 
simply be derived from the raw magnetic gradient data 
5 which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

Preferably the trend of the gradient of the ship bias is 
10 determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

Preferably, the further processing comprises the 
application of the corrected gradient data over sample 
15 intervals , integration of the corrected gradient data into 
total magnetic intensity data, and application of a low 
passed filter to the integrated total magnetic intensity 
data to provide the data output. 

2 0 Preferably the total magnetic intensity is obtained by 

integrating the corrected gradient data in accordance with 
the following equation: 

TMI= J) [GcAx,<0]+M, w , 

25 wherein G c is the corrected gradient data 

obtained after subtraction of the ship bias trend from the 
raw gradient data, and 

M Bt at is the total magnetic intensity value at a 
point representing the. start of the survey or at a 

3 0 location where a survey line and a tie line cross each 

other . 

The values of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter. 

35 

Preferably the gradient of the ship bias is calculated 
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according to the following equation: 

_ M,(x)-M,(x) [D(f,)-D(f 2 )] + [M ft (r,)-M fc (f 2 )] 

^bias — ™T 

A/ AZ 

5 where : 

MKx) = M,(x) + D(fi) + W,) 
MX*) = M«(x) + D(f 2) + M*(r 2 ) 

where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Me(x), diurnal variation D(tx), ship bias M b (tx) caused by 
ship-induced field, sensor drifting and sensor 
15 perturbation, all at time ti and along track distance x, 
and at some later time t 2 , Al is the distance between the 
sensors, D(t 2 ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
2 0 rear sensor M r . 

Xn one embodiment of the invention the first and second 
sensors which are towed behind the ship are included in a 
group of three or more towed sensors. 

25 

Xn this embodiment, preferably the number of sensors towed 
behind the ship comprises three sensors. 

According to this embodiment data from any two of the 
30 sensors is used to provide the raw magnetic gradient data. 

However, the embodiment may obtain raw gradient data from 
all three sensors and determine the trend of the gradient 
of the ship bias detected by all three sensors. 

35 

The invention may also be said to reside in a method of 
subsurface exploration to determine viability of drilling 
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in a marine environment: by considering magnetic data 
relating to the environment, and wherein the magnetic data 
has been obtained by: 

obtaining raw magnetic gradient data from the 

5 sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 
10 processing the corrected gradient data to provide 

a data output. 

Preferably, obtaining the magnetic data includes the step 
of determining the gradient of the ship bias from data 
15 obtained by the sensors, and determining the trend of the 
gradient of the ship bias from the gradient of the ship 
bias . 

Preferably the raw magnetic gradient data is calculated as 
2 0 the measured magnetic signal at the leading sensor less 
the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors* 

The estimate of the gradient of the ship bias may in fact 

2 5 simply be derived from the raw magnetic gradient data 

which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

30 Preferably the trend of the gradient of the ship bias is 
determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8* 

Preferably the processing comprises the application of the 

3 5 corrected gradient data over sample intervals, integration 

of the corrected gradient data into total magnetic 
intensity data, and application of a low passed filter to 
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the integrated total magnetic intensity data to provide 
the data output. 

Preferably the total magnetic intensity is obtained by 
5 integrating the corrected gradient data in accordance with 
the following equation: 

wherein G c is the corrected gradient data 
10 obtained after subtraction of the ship bias trend from the 
raw gradient data, and 

M 8 tat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
15 other. 

The values of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter. 

2 0 Preferably the gradient of the ship bias is calculated 
according to the following equation: 

Gbias = - = - 



25 where: 



30 



Mj{x) = M e {x)+D{t x )+M b (t x ) 
MM = M e (x) + D(/ 2 ) + M b (t 2 ) 



where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Mo(x) # diurnal variation D(ti), ship bias Mb(ti) caused by 
ship-induced field, sensor drifting and sensor 
3 5 perturbation, all at time t x and along track distance x, 
and at some later time t 2 , Al is the distance between the 
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sensors, D(t 2 ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
rear sensor M r . 

5 

In one embodiment of the invention the first and second 
sensors which are towed behind the ship are included in a 
group of three or more towed sensors* 

10 In this embodiment, preferably the number of sensors towed 
behind the ship comprises three sensors. 

According to this embodiment data from any two of the 
sensors is used to provide the raw magnetic gradient data. 

15 

However, the embodiment may obtain raw gradient data from 
all three sensors and determine the trend of the gradient 
of the ship bias detected by all three sensors. 

2 0 The invention may still further be said to reside in a 

method of drilling for a deposit in a marine environment, 
including: 

determining the location of drilling from data 
which has been obtained and which indicates the possible 
25 existence of the deposit; and 

which location is also determined by magnetic 
data which has been obtained by towing magnetic sensors 
behind a ship, the magnetic data being processed by: 

obtaining raw magnetic gradient data from the 

3 0 sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 
subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; 
35 and 

processing the corrected gradient data to provide 
a data output. 
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Preferably the method Includes the step of determining an 
estimate of the gradient of the ship bias from the raw 
magnetic gradient data obtained by the sensors, and 
5 determining the trend of the gradient of the ship bias 
from that estimate of the gradient of the ship bias. 

Preferably the raw magnetic gradient data is calculated as 
the measured magnetic signal at the leading sensor less 
10 the measured magnetic signal at the trailing sensor, 
divided by the distance between the sensors. 

The estimate of the gradient of the ship bias may in fact 
sixqply be derived from the raw magnetic gradient data 
15 which is obtained from the sensors and which may include 
data relating to the gradient of the ship bias as well as 
data relating to other magnetic effects. 

Preferably the trend of the gradient of the ship bias is 
20 determined by applying the gradient of the ship bias data 
to a B-Spline filter with smoothness 0.8. 

Preferably the processing of the corrected gradient data 
comprises the application of the corrected gradient data 
25 over sample intervals, integration of the corrected 

gradient data into total magnetic intensity data, and 
application of a low passed filter to the integrated total 
magnetic intensity data to provide the data output. 

3 0 Preferably the total magnetic intensity is obtained by 

integrating the corrected gradient data in accordance with 
the following equation: 

TMI=£ [G t Ax,(/)]+A/ iM , 

35 wherein G c is the corrected gradient data 

obtained after subtraction of the ship bias trend from the 
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raw gradient data, and 

Mstat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
5 other. 

The values of TMI so calculated are then smoothed to 
remove residual noise using a low pass filter. 

10 Preferably the gradient of the ship bias is calculated 
according to the following equation: 

_ M f (x)-M r (x) _ [D(f t ) - D(t 2 )] + [M b (t))-M h (t 2 )] 
Al Al 

15 where: 

M£x) = M e {x)+D{t x )+M b {t x ) 
MM = M e (x) + D(r 2 ) + M b (t 2 ) 

20 

where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Me(x), diurnal variation D(t x ), ship bias M b (ti) caused by 
ship- induced field, sensor drifting and sensor 
25 perturbation, all at time ti and along track distance x, 
and at some later time t 2 , Al is the distance between the 
sensors, D(t 2 ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
30 rear sensor M r . 

In one embodiment of the invention the first and second 
sensors which are towed behind the ship are included in a 
group of three or more towed sensors. 



35 



In this embodiment, preferably the number of sensors towed 
behind the ship comprises three sensors. 
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According to this embodiment data from any two of the 
sensors is used to provide the raw magnetic gradient data. 

5 However, the embodiment may obtain raw gradient data from 
all three sensors and determine the trend of the gradient 
of the ship bias detected by all three sensors. 

Brief Description of the Drawings 
10 A preferred embodiment of the invention will be described, 
by way of example, with reference to the accompanying 
drawings in which: 

Figure 1 is a view showing a marine magnetic data 

survey; 

15 Figure 2 is a flow chart showing the preferred 

method of processing the marine magnetic gradient data 
according to the preferred embodiment of the invention; 

Figure 2A is a graph showing the function of a 
low pass filter used in the preferred embodiment; 

2 0 Figure 3A and Figure 3B are graphs showing data 

obtained according to the survey and as processed 
according to the preferred embodiment of the invention; 

Figure 4 is a representation of a survey using 
data created according to a prior art technique; 
25 Figure 5 is a diagram similar to Figure 4, but 

using data obtained according to the preferred embodiment 
of the present invention; and 

Figure 6 is a view of a further embodiment of the 

invention. 

30 

Description of the Preferred Embodiment 

With reference to Figure 1, a marine survey to acquire 
magnetic data is illustrated. A ship 10 tows first 
magnetometer Ml and second magnetometer M2 by a tow line 

3 5 12. The ship travels along survey line 14 and the line 

16, which are dotted, represent the ship-induced magnetic 
field. 
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In Figure 1, it is assumed the ship is in the southern 
hemisphere and is heading north-northwest, and the sensors 
are drifting towards the northeast. Zn this case, the 
5 recorded gradient data between sensors Ml and M2 will have 
a larger distortion than if the sensors were drifting 
towards the southwest. 

Ideal conditions would assume that the sensors Ml and M2 
10 are towed directly behind the ship on the survey line and 
that therefore, after a given time period, the" sensor M2 
would occupy the same position as was previously occupied 
by the sensor Ml. However, due to drifting of the sensors 
because of ocean currents and the like, this ideal 
15 situation in general cannot be achieved. Nevertheless, 
the processing methods used to process magnetic data do 
assume that the sensor M2 will occupy the same position as 
the sensor Ml, previously occupied, but at a later time. 
The processing method according to the preferred 

2 0 embodiment enables ship bias to be accounted for much more 

accurately than any prior art techniques and therefore 
enables the sensors to be towed closer to the ship and to 
be separated by a smaller distance. This therefore 
provides much greater accuracy in the measurements because 
25 the amount of drift is not as great due to the shorter tow 
line and also because the separation between the sensors 
is smaller, the likelihood of the rear sensor occupying 
the same position as the front sensor at a later period in 
time is greatly increased. 

•30 

At any time ti and along- track distance x, the front sensor 
measures magnetic field - Mf consisting of the environment 
field Me(x), the diurnal variation D(ti), ship bias 2>&(ti) 
caused by ship induced field, sensor drifting and sensor 

3 5 perturbation. 



Mj{x) = M e (x) + Z>(/,) + Af*(/,) 
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At some later time t 2 , the rear sensor makes a measurement 
at the same spatial point: 

5 Mr(x) = M e (x) + D(t 2 ) + M b {t 2 ) 

The difference between the two measurements at the same 
location is: 

10 Af/x) - MM = [D(t x ) - D(t 2 )) + [M b {t x ) - M b Q 2 )) 

Ship bias gradient (Step 2, Figure 2) can be expressed as: 

= M f (x)~M r (x) ^ [D(r 1 ) - D(t 2 )] + [Af ^(r, ) - A/^ (r 2 )] 

15 

It should be mentioned that the calculated ship bias 
gradient is still affected by diurnal variation. In 
practice we notice that the sea current drifting is a kind 
of long wave-length variation so ship bias should be 
20 coincident with such current drifting. Step 2 in Figure 

2, which is the calculation of the gradient of the ship 
bias, is actually an estimate of the gradient of the ship 
bias which is obtained from the raw data collected by the 
sensors Ml and M2 . The data collected by the sensors will 

25 include raw magnetic gradient data which contains many 

signals including environmental magnetic signals, diurnal 
signals and the ship bias, as well as instrument bias and 
drift. The trend of the gradient of the ship bias is 
obtained from this estimate by applying the estimate of 

3 0 the gradient of the ship bias to a B-Spline filter, as 

will be described in more detail hereinafter. Therefore 
the trend of the gradient of the ship bias (Gtrend) (Step 

3, Figure 2) is used for the correction along a track 
line. 

35 



The gradient correction of the ship bias (G f )can be 
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expressed as follows: 

G C = G- Gtrend (Step 4, Figure 2) 

5 In this equation, 6 is the raw magnetic gradient data. 

As said before such drifting can cause large bias effects. 
Consequently, when the sensor drifting and perturbation 
take place simultaneously the bias zig-zag will be around 
10 its mean values. The variable mean value of the bias 

along a survey line is regarded as bias trend r since the 
high frequencies of sensor perturbation only create random 
noise around the bias trend and after integration its 
effect can be eliminated from raw gradient data. 

15 

After removal of ship bias effects TMI (Total Magnetic 
Intensity) can be calculated (Step 6) through integration 
of magnetic gradient data: 

N 

z 

20 TMI= " [GcAxi(t)]+Mstat 

Here, Axj.(t) is the sampling distance along the survey 
line. M s t*t is the TMI value at a point of survey start or 
at the location where the survey line and a tie line cross 

2 5 each other. 

The TMI value calculated in Step 6 is then smoothed using 
a low pass filter (Step 7) to remove any features which 
have a greater rate of change of TMI with distance than is 

3 0 expected in the particular survey area. An example of the 

action of this filter is provided in Figure 2A in which 
trace 50 is the smooth TMI curve and trace 51 is the TMI 
data prior to smoothing. 

3 5 The data output obtained at Step 7 may include line 

levelling and data griding to provide a final data output. 
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Figures 3A and 3B represent an actual example of the 
preferred embodiment of the invention which relates to raw 
data obtained from a known region. 

5 

In Figure 3A, the trace 20 represents the total ship bias. 
The line 21 represents the bias trend and it cam be seen 
from the left hand side of the trace 20 to the right hand 
side of the trace 20, the bias trend 21 changes non- 
10 uniformly about the value of -0.08 on the graph of Figure 
3A. A particular value for the bias trend can therefore 
be obtained for various intervals and this value can be 
subtracted from the raw data in order to provide the 
corrected data. The trend of the gradient of the ship 
15 bias is preferably determined by calculating the gradient 
of the bias data 20 and then applying that gradient data 
to a B-Spline filter with smoothness 0.8 in order to 
provide the representation of the ship bias 21. 

2 0 In the preferred embodiment of the invention, the gradient 

of the bias data 20 is calculated simply by using the raw 
magnetic gradient data obtained from the sensors because 
when that data is applied to the B-Spline filter, all that 
remains is the gradient of the ship bias component. The 
25 ship bias 21 can then be subtracted from the calculated 
raw gradient. 

The trend of the gradient of the ship bias is a non-linear 
function which is represented by the line 21 in Figure 3A. 

3 0 The line is a measure of how the bias changes over time as 

the sensors are towed behind the ship* As is clearly seen 
from line 21, the trend is not constant or merely an 
average, but rather represents the fluctuation of the bias 
gradient and at some times, is above the value of 0.08 and 
35 sometimes below that value. As noted above, the trend is 
determined by applying the gradient of the ship bias to a 
B-Spline filter with smoothness 0.8. However, in other 
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embodiments, the filter could have different smoothness 
depending on the region from which data is collected and 
the nature of the data which is collected. Xn general, 
the function of the filter is to smooth the trace 20 so 
5 that some meaningful value of the bias trend at particular 
time periods can be obtained. Thus, the filter is 
effectively determining the peaks and troughs of the trace 
20 and fitting a curved line between those peaks and 
troughs, which gives a measure of how the gradient of the 
10 ship bias is changing with time as the sensors are towed 
behind the ship. 

In the example of data obtained from the known region, 
dashed line 30 represents total magnetic intensity data 

15 which is provided according to the prior art processing 
technicjue. The line 32 represents the integrated total 
magnetic intensity data from the bias corrected gradient 
data according to the preferred embodiment of the 
invention. The line 34 represents the diurnal variation 

20 at a station some 500 km away from the survey area and 

line 36 is the observed field data including the diurnal 
effect • 

It can be seen from the conventional processing technique, 
25 which gives line 30, false anomalies are provided which do 
not show up on the line 32, which is produced according to 
the present invention. 

Xn general, the results of the magnetic survey are 
.3,0 produced in a colour chart. Black and white 

representations of a colour chart are shown as Figure 4 in 
relation to conventional processing of data obtained from 
a known region and Figure 5 by a processing method 
according to the preferred embodiment of this invention. 

35 

The characteristics of the Igneous basement of the known 
region are well known and it can be seen that data 
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processed according to the present invention produces a 
much, clearer indication of actual magnetic structures than 
prior art technique which include considerable 
interference and data which may mislead analysers into 
5 considering that magnetic structures exist which are not 
actually there, or that magnetic structures do not exist 
when a magnetic structure actually is present. 

Thus, the preferred embodiment of the present invention 
10 provides data which more accurately reflects the 

likelihood of magnetic structures which can tlien be used 
as a basis for determining the viability of a drilling 
operation in an exploration or recovery program. 

15 Figure 6 shows a second embodiment of the invention in 

which like reference numerals indicate like components to 
those described with reference to Figure 1. In this 
embodiment, three magnetometers M 1# M 2 and M 3 are towed by 
the ship 10 on tow line 12. This embodiment therefore 

20 builds some redundancy into the system in that, should one 
of the sensors fails, two sensors are still available in 
order to provide the required gradient data, and therefore 
there is considerably less likelihood that a survey will 
be wasted should a ship complete the survey only to find 

25 that one of the magnetometers has not operated properly. 
In this embodiment, preferably the distance between the 
magnetometers Mi and M 2 , and the distance between the 
magnetometers M 2 and M 3 is about 15 metres. The towing 
distance between the ship 10 and the first magnetometer M x 

30 is preferably about 150 metres or less. 

This embodiment provides the further advantage that any 
group of two sensors can be used to provide the gradient 
data, and therefore the gradient data could be provided by 
3 5 the magnetometers M x and M 2 , the magnetometers M 2 and M 3 , or 
the magnetometers Mi and M 3 . This embodiment provides the 
further advantage that all three magnetometers could be 
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used to provide the data to enable the calculation of the 
raw gradient and the gradient of the ship bias. The use 
of three magnetometers may improve accuracy because of the 
increased amount of data which is collected. 

5 

Obviously, more than three magnetometers could be used if 
desired, thereby increasing the amount of data collected 
and reducing the likelihood that a survey will be wasted 
because at least two of the magnetometers are not operated 
10 properly. 



Since modifications within the spirit and scope of the 
invention may readily be effected by persons skilled 
within the art, it is to be understood that this invention 
15 is not limited to the particular embodiment described by 
way of example hereinabove. 
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Claims 



!• A method of processing marine magnetic data 

obtained by towing first and second spaced apart sensors 
5 behind a ship, the method including: 

obtaining raw magnetic gradient data from the 

sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 
10 subtracting the trend from the raw magnetic 

gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
a data output. 

15 2. The method of claim 1 including the step of 

determining an estimate of the gradient of the ship bias 
from the raw magnetic gradient data obtained by the 
sensors, and determining the trend of the gradient of the 
ship bias from that estimate of the gradient of the ship 

2 0 bias. 



3. The method of claim 1 wherein the raw magnetic 

gradient data is calculated as the measured magnetic 
signal at the leading sensor less the measured magnetic 
25 signal at the trailing sensor, divided by the distance 
between the sensors. 



4. The method of claim 1 wherein the trend of the 

gradient of the ship bias is determined by applying the 
3 0 gradient of the ship bias data to a B-Spline filter with 
smoothness 0.8. 



5. The method of claim 1 wherein the processing 

comprises the application of the corrected gradient data 
3 5 over sample intervals, integration of the corrected 

gradient data into total magnetic intensity data, and 
application of a low pass filter to the integrated total 
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magnetic intensity data to provide the data output* 



6. The method of claim 5 wherein the total magnetic 
intensity is obtained by integrating the corrected 

5 gradient data in accordance with the following equation: 

TMI=£ [GcAxm+Msw 

M 

wherein G c is the corrected gradient data 
obtained after subtraction of the trend of the gradient of 
10 the ship bias from the raw gradient data, and 

Mstat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
other . 

15 

7. The method of claim 1 wherein the gradient of the 
ship bias is calculated according to the following 
equation: 

20 G M f {x)-M r {x) ^ [D(t x )-D{t 2 )] + [M h {t x )-M b {t 2 y\ 

bias Al ^ Al 

where : 



25 



Af/^) = M tf (x) + D(r,) + ^(f|) 
MM = MM + D{t 2 ) + M b (t 2 ) 



where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
30 Me(x), diurnal variation D(ti), ship bias M b (ti) caused by 
ship-induced field, sensor drifting and sensor 
perturbation, all at time t x and along track distance x, 
and at some later time t 2 # Al is the distance between the 
sensors, D(t 2 ) is the diurnal variations sensed by the rear 
3 5 sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
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rear sensor M r . 

8. The method, of claim 1 wherein the first and 
second sensors which are towed behind the ship are 

5 included in a group of three or more towed sensors. 

9. The method of claim 8 wherein the number of 
sensors towed behind the ship comprises three sensors* 

10 10. A method of obtaining gradient data for an 

exploration program, the method including: 

towing first and second sensors behind a ship 
along predetermined survey lines; 

obtaining raw magnetic gradient data from the 

15 sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 
2 0 processing the corrected gradient data to provide 

a data output* 

11. The method of claim 10 wherein the method 

includes the step of determining an estimate of the 

2 5 gradient of the ship bias from the raw magnetic gradient 

data obtained by the sensors, and determining the trend of 
the gradient of the ship bias from that estimate of the 
gradient of the ship bias. 

3 0 12. The method of claim 10 wherein the raw magnetic 

gradient data is calculated as the measured magnetic 
signal at the leading sensor less the measured magnetic 
signal at the trailing sensor, divided by the distance 
between the sensors. 

35 

13. The method of claim 10 wherein the trend of the 

gradient of the ship bias is determined by applying the 
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gradient of the ship bias data to a B-Spline filter with 
smoothness 0.8. 



14. The method of claim 10 wherein the further 

5 processing comprises the application of the corrected 
gradient data over sample intervals, integration of the 
corrected gradient data into total magnetic intensity 
data, and application of a low passed filter to the 
integrated total magnetic intensity data to provide the 
10 data output. 

15. The method of claim 14 wherein the total magnetic 
intensity is obtained by integrating the corrected 
gradient data in accordance with the following equation: 

15 

TMI=X [GcAx,</)]+M JftJ/ 

wherein G c is the corrected gradient data 
obtained after subtraction of the ship bias trend from the 
raw gradient data, and 
20 M 8t at is the total magnetic intensity value at a 

point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
other . 

25 16. The method of claim 10 wherein the gradient of 

the ship bias is calculated according to the following 
equation: 



M,(x)-M,(x) [Dit^-Djt^ + WM-M^)) 
A/ A/ 



30 

where : 



35 



Mjix) = M e (.x) + D(/,) + M„(t0 
MM = M e (x) + D(r 2 ) + M b (t 2 ) 



WO 03/052460 



PCT/AU02/01669 



- 27 - 

where M f is the magnetic field as measured by the 
front sensor and which consists of the environmental field 
Ma(x), diurnal variation D(ti)/ ship bias M b (t!) caused by 
ship-induced field, sensor drifting and sensor 
5 perturbation, all at time t x and along track distance x, 
and at some later time t 2 , Al is the distance between the 
sensors, D(t 2 ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
10 rear sensor M r . 



17. The method of claim 10 wherein the first and 
second sensors which are towed behind the ship are 
included in a group of three or more towed sensors. 

15 

18. The method of claim 17 wherein the number of 
sensors towed behind the ship comprises three sensors. 



19. A method of subsurface exploration to determine 

20 viability of drilling in a marine environment by 

considering magnetic data relating to the environment, and 
wherein the magnetic data has been obtained by: 

obtaining raw magnetic gradient data from a 
plurality of sensors; 
25 determining the trend of the gradient of the ship 

bias detected by the sensors; 

subtracting the trend from the raw magnetic 
gradient data to obtain corrected gradient data; and 

processing the corrected gradient data to provide 
30 a data output. 



20. The method of claim 19 wherein the magnetic data 

includes the step of determining the gradient of the ship 
bias from data obtained by the sensors, and determining 
3 5 the trend of the gradient of the ship bias from the 
gradient of the ship bias. 
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21. The method of claim 19 wherein the raw magnetic 
gradient data is calculated as the measured magnetic 
signal at the leading sensor less the measured magnetic 
signal at the trailing sensor, divided by the distance 

5 between the sensors. 

22. The method of claim 19 wherein the trend of the 
gradient of the ship bias is determined by applying the 
gradient of the ship bias data to a B-Spline filter with 

10 smoothness 0.8. 

23. The method of claim 19 wherein the processing 
comprises the application of the corrected gradient data 
over sample intervals, integration of the corrected 

15 gradient data into total magnetic intensity data, and 

application of a low passed filter to the integrated total 
magnetic intensity data to provide the data output. 

24. The method of claim 23 wherein the total magnetic 
2 0 intensity is obtained by integrating the corrected 

gradient data in accordance with the following equation: 

TMI= £ [G c *Xi{t)]+M S!at 
i=i 

wherein G c is the corrected gradient data 
25 obtained after subtraction of the ship bias trend from the 
raw gradient data, and 

Mfltat is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
30 other. 

25. The method of claim 19 wherein the gradient of 
the ship bias is calculated according to the following 
equation: 
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M,(x)-M,(x) _ ) - D(t 3 )] + [M„ (f , ) - M„ (t t )] 
Gwas - 

where : 

5 Mj(x) = M e (x) + D(ri) + M 6 (ri) 

M*x) = (x) + D(/ 2 ) + Afcfe) 

where M f is the magnetic field as measured by the 
10 front sensor and which consists of the environmental field 
Me(x), diurnal variation DCtx), ship bias M b (tiT caused by 
ship-induced field, sensor drifting and sensor 
perturbation, all at time ti and along track distance x, 
and at some later time t 2 , Al is the distance between the 
15 sensors, D(t 2 ) is the diurnal variations sensed by the rear 
sensor, and 

M b (t 2 ) is the ship bias at time t 2 detected by the 
rear sensor M r . 



2 0 26. The method of claim 19 wherein the number of 

sensoirs towed behind the ship comprises three sensors. 

27. A method of drilling for a deposit in a marine 

environment , including : 
25 determining the location of drilling from data 

which has been obtained and which indicates the possible 
existence of the deposit; and 

which location is also determined by magnetic 
data which has been obtained by towing magnetic sensors 

3 0 behind a ship, the magnetic data being processed by: 

obtaining raw magnetic gradient data from the 
sensors; 

determining the trend of the gradient of the ship 
bias detected by the sensors; 
3 5 subtracting the trend from the raw magnetic 

gradient data to obtain corrected gradient data; 
and 
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processing the corrected gradient data to provide 
a data output. 

28. The method of claim 27 wherein the method 
5 includes the step of determining an estimate of the 

gradient of the ship bias from the raw magnetic gradient 
data obtained by the sensors, and determining the trend of 
the gradient of the ship bias from that estimate of the 
gradient of the ship bias. 

10 

29. The method of claim 27 wherein the xfaw magnetic 
gradient data is calculated as the measured magnetic 
signal at the leading sensor less the measured magnetic 
signal at the trailing sensor, divided by the distance 

15 between the sensors. 

30. The method of claim 27 wherein the trend of the 
gradient of the ship bias is determined by applying the 
gradient of the ship bias data to a B-Spline filter with 

2 0 smoothness 0.8. 

31. The method of claim 27 wherein the processing of 
the corrected gradient data comprises the application of 
the corrected gradient data over sample intervals, 

25 integration of the corrected gradient data into total 

magnetic intensity data, and application of a low passed 
filter to the integrated total magnetic intensity data to 
provide the data output. 

30 32. The method of claim 31 wherein the total magnetic 

intensity is obtained by integrating the corrected 
gradient data in accordance with the following equation: 

TMI= £ [G<Ax,<0]+M, w , 

/a| 

35 wherein G c is the corrected gradient data 

obtained after subtraction of the ship bias trend from the 
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raw gradient data, and 

M st at is the total magnetic intensity value at a 
point representing the start of the survey or at a 
location where a survey line and a tie line cross each 
5 other* 

33. The method of claim 27 wherein the gradient of 

the ship bias is calculated according to the following 
equation: 

10 

_ M f (x)-M,(x) _ [0(r, ) - D(r 2 )] + [M b (tf)-M h (t 2 )] 
where : 

1 5 Mj(x) = M e (x) + D(r0 + M b (t x ) 

MAx) = M e (x) + D(t 2 ) + M b (t 2 ) 

where M f is the magnetic field as measured by the 
20 front sensor and which consists of the environmental field 
Mo(x), diurnal variation Dfti), ship bias M b (t 2 ) caused by 
ship- induced field, sensor drifting and sensor 
perturbation, all at time t x and along track distance x, 
and at some later time t 2 , Al is the distance between the 

2 5 sensors, D(t 2 ) is the diurnal variations sensed by the rear 

sensor, and 

M b (t 2 ) is the ship bias at time t 3 detected by the 
rear sensor M r . 

3 0 34. The method of claim 27 wherein the first and 

second sensors which are towed behind the ship are 
included in a group of three or more towed sensors. 

35. The method of claim 34 wherein the number of 

35 sensors towed behind the ship comprises three sensors. 



36. 



The method of claim 35 wherein data from any two 
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of the sensors is used to provide the raw magnetic 
gradient data. 
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Calculate raw gradient — Step 1 



Caculate gradient of ship bias — Step 2 



I 



Calculate trend of the gradient of the ship bias — Step 3 
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Subtract bias trend from raw gradient — Step 4 



I 



Apply corrected gradient «. 5 

over sample intervals °' ' 



1 



Integration of gradient into TMI — Step 6 
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Apply low-pass filter to the 
integrated TMI 



— Step 7 
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Data output 
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